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^ ■ Abstract 

■ This paper is concerned with the existence of the simultaneous route- and- departure 

(-H , choice dynamic user equilibrium (SRDC-DUE) in continuous time. The SRDC-DUE prob- 

lem was formulated as an infinite-dimensional variational inequality in Friesz et al. (1993). 
In deriving our existence result, we employ the generalized Vickrey model (GVM) intro- 
duced in Han et al. (2012a, 2012b) to formulate the underlying network loading problem. 
As we explain, the GVM corresponds to a path delay operator that is provably strongly 
continuous on the Hilbert space of interest. Finally, we provide the desired SRDC-DUE 
. existence result for general constraints relating path flows to a table of fixed trip volumes 

without invocation of a priori bounds on the path flows. 
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1 Introduction 



I In this paper we shall consider dynamic traffic assignment (DTA) to be the positive (de- 

scriptive) modeling of time-varying flows of automobiles on road networks consistent with 
established traffic flow theory and travel demand theory. Dynamic User Equilibrium (DUE) 
, is one type of DTA wherein effective unit travel delay for the same purpose is identical for all 

^ ■ utilized path and departure time pairs. The relevant notion of travel delay is effective unit 



travel delay, which is the sum of arrival penalties and actual travel time. For our purposes 
in this paper, DUE is modeled for the within-day time scale based on fixed travel demands 
established on a day-to-day time scale. 

In the last two decades there have been many efforts to develop a theoretically sound 
formulation of dynamic network user equilibrium that is also a canonical form acceptable to 
scholars and practitioners alike. DUE models tend to be comprised of four essential sub- 
models: 

1. a model of path delay; 

2. flow dynamics; 
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3. flow propagation constraints; 

4. a path/departure-time choice model. 

Furthermore, analytical DUE models tend to be of two varieties: (1) route choice (RC) user 
equilibrium (Friesz et al. 1989, Merchant and Nemhauser 1978a, 1978b, Mounce 2006, Smith 
and Wisten 1995, Zhu and Marcotte 2000) ; and (2) simultaneous route- and- departure choice 
(SRDC) dynamic user equilibrium (Friesz et al. 1993, 2010, 2012, Ran et al. 1996, Wie 
et al. 2002). For both types of DUE models, the existence of a dynamic user equilibrium 
in continuous time remains a fundamental issue. A proof of DUE existence is a necessary 
foundation for qualitative analysis and computational studies. In this paper, we provide 
a DUE existence result for the SRDC DUE problem when it is formulated as an infinite- 
dimensional variational inequality of the type presented in Friesz et al. (1993). In this paper, 
in order to establish a DUE existence result, we study the network loading problem based 
on the generalized Vickrey model (GVM) proposed in Han et al. (2012a, 2012b). All of our 
results presented in this paper are more general than any obtained previously for DUE when 
some version of the point queue model is employed. 

1.1 Formulation of the SRDC user equilibrium 

There are two essential components within the RC or SRDC notions of DUE: (i) the mathe- 
matical expression of Nash-like equilibrium conditions, and (ii) a network performance model, 
which is, in effect, an embedded network loading problem. The embedded network loading 
problem captures the relationships among arc entry flow, arc exit flow, arc delay and path 
delay for any path departure rate trajectory. Note that, by studying the embedded network 
loading problem based on the GVM, we are not suggesting or employing a sequential so- 
lution paradigm for DUE. If the reader believes a sequential perspective is implicit in our 
work, he/she has not fully understood the mathematical presentation provided in this and 
subsequent sections. 

There are multiple means of expressing the Nash-like notion of a dynamic equilibrium, 
including the following: 

1. a variational inequality (Friesz et al. 1993; Smith and Wisten 1994, 1995) 

2. an equilibrium point of an evolution equation in an appropriate function space (Mounce 
2006; Smith and Wisten 1995) 

3. a nonlinear complementarity problem (Wie et al. 2002; Han et al. 2011) 

4. a differential variational inequality (Friesz et al. 2001, 2010, 2012; Friesz and Mookherjee 
2006); and 

5. a differential complementarity system (Pang et al. 2011). 

The variational inequality representation is presently the primary mathematical form employed 
for both RC and SRDC DUE. The most obvious approach to establishing existence for any of 
the mathematical representations mentioned above is to convert the problem to an equivalent 
fixed point problem and then apply Brouwer's fixed point existence theorem. Alternatively, 
one may use an existence theorem for the particular mathematical representation selected; it 
should be noted that most such theorems are derived by using Brouwer's famous theorem. So, 
in effect, all proofs of DUE existence employ Brouwer's fixed point theorem, either implicitly 
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or explicitly. One statement of Brouwer's theorem appears as Theorem 2 of Browder (1968). 
Approaches based on Brouwer's theorem require the set of feasible path flows (departure rates) 
under consideration to be compact and convex in a Banach space, and typically involve an a 
priori bound on each path flow. 

We also wish to point out that this paper employs much more general constraints relating 
path flows to a table of fixed trip volumes than has been previously considered when studying 
SRDC DUE. Moreover, in our study of existence, we do not invoke a priori bounds on the 
path flows to assure boundedness needed for application of Brouwer's theorem. That is, a 
goal of this paper is to investigate the existence of DUE without making the assumption of a 
priori bounds for departure rates. Note should be taken of the following fact: the boundedness 
assumption is less of an issue for the RC DUE by virtue of problem formulation; that is, for 
RC DUE, the travel demand constraints are of the following form: 

J2 hp{t) = R,j{t), Vt, y{i,j)eW (1.1) 

where W is the set of origin-destination pairs, Vij is the set of paths connecting G W 
and hp{t) is the departure rate along path p. Furthermore, Rijit) represents the rate (not 
volume) at which travelers leave origin i with the intent of reaching destination j at time t; 
each such trip rate is assumed to be bounded from above. Since (jl.ip is imposed pointwise and 
every path flow hp is nonnegative, we are assured that each h = {hp : G yV,p G Vij) are 
automatically uniformly bounded. On the other hand, the SRDC user equilibrium imposes 
the following constraints on path flows: 

r Kit)dt = Qij, V(i, j)GW (1.2) 

where Qij G is the volume (not rate) of travelers departing node i with the intent of 
reaching node j. The integrals in (jl.2p are interpreted as Legesgue; hence, (jl.2p alone is not 
enough to assure bounded path flows. This observation has been the major hurdle to proving 
existence without the a priori invocation of bounds on path flows. In this paper, we will 
overcome this difficulty through careful analysis of the GVM and by investigating the effect 
of user behavior in shaping network flows, in a mathematically intuitive yet rigorous way. 

1.2 Importance of the path delay operator 

Clearly another key component of continuous-time DUE is the path delay operator, typically 
obtained from dynamic network loading (DNL), which is a subproblem of a complete DUE 
mode0. Any DNL must be consistent with the established path flows and link delay model, 
and DNL is usually performed under the fir st-in- first- out (FIFO) rule. The properties of 
the delay operator are critical to proving existence of a solution to the infinite-dimensional 
variational inequality used to express DUE. In Zhu and Marcotte (2000), using the link delay 
model introduced by Friesz et al. (1993), the authors showed weak continuity of the path delay 

^Note that, by refering to the network loading procedure, we are neither employing nor suggesting a se- 
quential approach to the study and computation of DUE. Rather a subset of the equations and inequalites 
comprising a complete DUE model may be grouped in a way that identifies a traffic assignment subproblem 
and a network loading subproblem. Such a grouping and choice of names is merely a matter of convenient 
language that avoids repetitive reference to the same mathematical expressions. Use of such language does not 
alter the need to solve both the assignment and loading problems consistently and, thus, simultaneously. A 
careful reading of the mathematical presentation made in subsequent sections makes this quite clear. 
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operator under the assumption that the path flows are a priori bounded. Then' continuity 
result is superceded by a more general result proven in this paper: the path delay operator 
of interest is strongly continuous without the assumption of boundedness. Strong continuity 
without boundedness is central to our proof of existence in the present paper. 

In this paper, as a foundation for DNL, we will consider Vickrey's model of congestion 
first introduced in Vickrey (1969) and later studied by Han et al. (2012a, 2012b). Vickrey's 
model for one link is primarily described by an ordinary differential equation (ODE) with 
discontinuous right hand side. Such irregularity has made it difficult to analyze Vickery's 
model in continuous time. Fortunately, in this paper, we will be able to take advantage of the 
reformulation proposed in Han et al. (2012a, 2012b), then prove the strong continuity of the 
path delay operator without boundedness of the path flows. This will provide a quite general 
existence proof for SRDC-DUE based on the GVM. 

1.3 Organization 

The balance of this paper is organized as follows. Section [2] provides essential mathematical 
background on the concepts that will be used in the paper. Section [3] briefly reviews the formal 
definition of dynamic user equilibrium and its formulation as a variational inequality. Section 
m recaps the generalized Vickrey model (GVM) originally put forward by Han et al. (2012a, 
2012b). Section [5] formally discusses the properties of the effective delay operator. The main 
result of this paper, the existence of an SRDC-DUE when the GVM informs network loading 
is established in Theorem 15.71 of Section 15.31 

2 Mathematical preliminaries 

A topological vector space is one of the basic structures investigated in functional analysis. 
Such a space blends a topological structure with the algebraic concept of a vector space. The 
following is a precise definition: 

Definition 2.1. (Topological vector space) A topological vector space X is a vector space 
over a topological field ¥ (usually the real or complex numbers with their standard topologies) 
which is endowed with a topology such that vector addition X x X ^ X and scalar multipli- 
cation F X X — 7- X are continuous functions. 

As a consequence of Definition 12.11 normed vector spaces, and therefore all Banach spaces 
and Hilbert spaces, are examples of topological vector spaces. Also important is the notion of 
a seminorm: 

Definition 2.2. (Seminorm) A seminorm on a vector space X is a real-valued function p 
on X such that 

(a) p{x + y) < p{x) +p{y) 

(b) p{a x) = \a\ p{x) 

for all X and y in X and all scalars a. 

Definition 2.3. (Locally convex space) A locally convex space is defined to be a vector 
space X along with a family of seminorms {pi}iex on X. 

As part of our review we make note of the following essential knowledge: 
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Fact 1. The space of square-integrable real- valued functions on a compact interval [a, b], 
denoted by /2^([a, b]), is a locally convex topological vector space. 

Fact 2. The m-fold product of the spaces of square-integrable functions C^{[a, 6])"^ is a 
locally convex topological vector space. 

Definition 2.4. (Dual space) The dual space X* of a vector space X is the space of all 

continuous linear functions on X 

Another key property we consider without proof is: 

Fact 3. The dual space of L^ip,) for 1 < p < oo has a natural isomorphism with L''{p,) where 
q is such that 1/p + l/q = 1. In particular, the dual space of C^{[a, b]) is again £^([a, b]). 

The key foundation for analysis of existence is the following theorem given in Browder (1968): 

Theorem 2.5. Let K be a compact convex subset of the locally convex topological vector space 
E, T a continuous (single-valued) mapping of K into E* . Then there exits uq in K such that 



VI{F, U) (2.3) 



T{uo), UQ-u) > 
for all u K . 

Let us now give the formal definition of a variational inequality in a topological setting: 

Definition 2.6. (Infinite- Dimensional Variational inequality) Let V be a topological 
vector space and F : U x 3?^. — )■ V, where U C V. The infinite- dimensional variational 
inequality is posed as the following problem 

find u* eU 
{F{u),u-u*) > y ueU 

Definition 2.7. ( Compactness of subspaces) A subset K of a topological space X is called 
compact if for every arbitrary collection {C/a}^g^ of open subsets of X such that 

K C (jUa 

there is a finite subset I of A such that 

K C \JUi 

Definition 2.8. (Sequential compactness) A topological space is sequentially compact if 
every sequence has a convergent subsequence. 

An outgrowth of the concepts and results given above, the following fact is stated without 
proof: 

Fact 4. In metric space (hence topological vector space), the notions of compactness and 
sequential compactness are equivalent. 

The final bit of specialized knowledge about topological vector spaces that we shall need is 
the following: 

Definition 2.9. (Weak convergence in Hilbert space) A sequence of points {xn} in a 
Hubert space 7i is said to be convergent weakly to a point x € H, denoted as Xn ^ y if 

{xn, v) (x, y) 

for all y GTi, where (•, •) is the inner product on the Hilbert space. 



5 



3 Continuous-time dynamic user equilibrium 



In this section, we will assume the time interval of interest is 

[to, tf]c^\ 

The most crucial component of the DUE model is the path delay operator, which provides 
the time to traverse any path p per unit of flow departing from the origin of that path. The 
delay operator is denoted by 

Dp{t, h) y p£V 

where V is the set of all paths employed by network users, t denotes the departure time, and 
h is a vector of departure rates. Throughout the rest of the paper, we stipulate that 

where ^>C^[to, tf]j denotes the non-negative cone of the iPj-fold product of the Hilbert 
space C'^[to, tj] of square-integrable functions on the compact interval [to, t/]. The inner 
product of the Hilbert space ( [to , t/] ) is defined as 



{u{s)) v{s)ds (3.4) 



to 



where the superscript T denotes transpose of vectors. Moreover, the norm 

||ti||^2 = {u, u)^^"^ (3.5) 

is induced by the inner product ()3.4p . 

Next, we need to consider a more general notion of travel cost that will motivate on-time 
arrivals. To this end, for each p (z V, we introduce the effective unit path delay operator 

: [to, tf] X (£^([to, tf])] — )• and define it as follows: 



^p{t, h) = Dp{t, h) + F\t + Dp{t, h) - TaJ (3.6) 

where J-'(-) is the penalty for early or later arrival relativerelative to the target arrival time 
Ta- Note that, for convenience, Ta is assumed to be independent of destination. However, 
that assumption is easy to relax, and the consequent generalization of our model is a trivial 
extension. We interpret ^p(t, h) as the perceived travel cost of driver starting at time t on 
path p under travel conditions h. Presently, our only assumption on such costs is that for each 

h G (£^([to, tf])j , the vector function ^(•, h) : [to, t/] ^^++ is measurable and stictly 



positive. The assumption of measurability was used for a measure theory-based argument in 
Friesz et al. (1993). Later in this paper, we shall discuss other properties of this operator, 
such as continuity on a Hilbert space. The continuity of effective delay is crucial for applying 
the general theorems in Browder (1968), especially Theorem 12.51 stated above. 

To support the development of a dynamic network user equilibrium model, we introduce 
some additional constraints. Foremost among these are the flow conservation constraints 

^ r hp{t) dt = Qij V (i, j) G W (3.7) 

J to 



p&p, 
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where Vij is the set of all paths that connect origin-destination (0-D) pair (i, j) G W, while 
W is the set of ah 0-D pairs. In addition, Qij is the fixed travel demand for 0-D pair (i, j). 
Using the notation and concepts we have thus far introduced, the set of feasible solutions for 
DUE when the effective delay operator ^(•, •) is given is 




Using a presentation very similar to the above, the notion of a dynamic user equilibrium 
in continuous time was first introduced by Priesz et al. (1993), who employ a definition 
tantamount to the following: 

Definition 3.1. (Dynamic user equilibrium). A vector of departure rates (path flows) 
h* ^ K is a dynamic user equilibrium if 

h; (t) > 0,p G V^J =^ [t, h* {t)] = Vij G V(i, j) G W (3.9) 

We denote the dynamic user equilibrium defined this way by DUE , A, [tQ,tf]). 

In the analysis to follow, we focus on the following infinite-dimensional variational inequality 
formulation of the DUE problem reported in Theorem 2 of Friesz et al. (1993). 

find /i* G A such that 

E r ^p{t,h*){hp-h;)dt>o \viU,A,[to,tf]) (3.10) 

pevJto ^ ' 

yh£A ) 

The variational inequality formulation Vl["i>, A, [to, tf]) expressed above subsumes almost all 
DUE models regardless of the arc dynamics or network loading models employed provided 

^p{t,h*) 

4 The dynamic network loading 

A key ingredient of the variational inequality formulation of the DUE (I3.10p is the effective 
delay operator ^'(t, •), which maps a vector of admissible departure rates to the vector of 
strictly travel costs associated with each route-and-departure-time choice. The problem of 
predicting time-varying network flows consistent with known travel demands and departure 
rates (path flows) is usually referred to as the dynamic network loading (DNL) problem. 
Since effective path delays are constructed from arc delays that depend on arc activity and 
performance, DNL is intertwined with the determination of effective delay operators. 

In this section we present a continuous-time DNL model. This model is based on a re- 
formulation of Vickrey's model (Vickrey 1969) that we refer to as the generalized Vickrey 
model (GVM); it was apparently first proposed in Han et al. (2012a, 2012b). The generalized 
Vickrey model determines arc exit flow and the arc traversal time from arc entry flow. This 
formulation not only leads to a simple and explicit computational scheme, but also makes it 
easier to conduct rigorous analyses of the arc delay operator and, hence, of the effective path 
delay operator 4'(t, •). 
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4.1 The generalized Vickrey model 

First introduced in Vickrey (1968), the Vickrey's model is based on two key assumptions: 
(i) vehicles have negligible sizes, and, therefore, any queue is of negligible size; and (ii) link 
travel time consists of a fixed travel time plus a congestion-related arc-traversal delay. Let us 
introduce the following notation: 



u{t) 
M 
Q{t) 

wit) 
T 

D{t) 



the flow profile into the link 

the flow capacity of the link 

the queue size 

the link exit flow 

the fixed free flow travel time 

the link traverse time of drivers entering the link at t 



Then the model is described by the following set of equations. 

fminHt-T),M} ,(t) = 

^ ^ \m q{t) / ^ ' 

^ = u{t-T)-w{t) (4.12) 

D{t)=T+^-^^ (4.13) 

Notice that (j4.1ip and (j4.12p amount to an ordinary differential equation (ODE) with a right 
hand side that is discontinuous in the state variable q{-). Such ODEs have been the main 
hurdle to further analysis and computation of this model in continuous time. In Han et 
al. (2012a, 2012b), a reformulation of Vickrey's model as a Hamilton- Jacobi equation was 
proposed and solved with a version of the Lax-Hopf formula (see Evans 2010 for more detailed 
on general Hamilton- Jacobi equation and Lax-Hopf formula). As a result, the solution to 
(j4.1ip - (|4.13|) was presented in closed form. Due to the space limitation, we will omit details 
of the derivation of explicit solution to Vickrey's model and refer the reader to Han et al. 
(2012a, 2012b) for details. 

Let us next introduce the cumulative entering vehicle count [/(•) and the exiting vehicle 
count W(-) at the entrance and exit of the link of interest, respectively. Furthermore, U{-) is 
assumed to be non-decreasing and left continuous. Notice that these latter assumption imply 
that the link entry flows can be unbounded and possibly contain the dirac delta function. 
In contrast, the Vickrey's original model requires that the entry flow to be at least Lebesgue 
integrable. As such, the GVM is more general than Vickrey's model. 

Using the notation introduced previously, an equivalent statement of (|4.1ip and (|4.12p is 
the following: 

W{t) = mm {U{t)-Mt} + M{t- to) (4.14) 

T<t — T 

q{t) = U{t)-Mt- min {U{t)-Mt} (4.15) 

T<t — T 

D{t) = T+'-^^ (4.16) 
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Note that all the unknowns of system (j4.14p - (j4.16|) may be explicitly stated in terms of the 
cumulative vehicle count U{-). The system (j4.14p - ()4.16p will serve as the mathematical for- 
mulation of link dynamics in the DNL, as we shall explain shortly. The system (j4.14p - (|4.16p 
may also be used for deriving mathematical properties of the effective path delay operator, as 
is demonstrated in Section [5Tl 



4.2 The network model 

It is straightforward to extend the generalized Vickrey model to a network, which is represented 
as a directed graph G{N, A), where and A are the set of nodes and arcs, respectively. In 
order to proceed, we introduce some additional notations. In particular, for each node v (z N, 
let I"" be the set of incoming links, the set of outgoing links. For each arc a € A, let Ua{t), 
Wa{t) be the entry flow and exit flow, respectively. The arc entry/exit flows are the sum of 
entry/exit flows associated with individual paths using this arc; that is 



u, 



a{t) = 5^<^ap<(i), Wait) = 5^<^ap«^SW, VoG^ (4.17) 



where 



^ap 



1 if arc a belongs to path p 

otherwise 



Let us also define the cumulative entering vehicle count Ua{t) and cumulative exiting vehicle 
count Wait), for each arc a. Similarly, each one is disaggregated into quantities associated 
with each path that uses this arc: 

Ua{t) = Y.6apUP{t), Wait) = ^^dapW^it), VcGA (4.18) 

p€V p€V 

The arc delay Dait) is the time taken to traverse the arc a when the time of entry is t. 
The arc exit time function Tait) is defined as ra(t) = t + Dait), that is, Tait) represents the 
exit time on arc a when the time of entry is t. 

For each group of drivers using the arc, the ratio of their arrival and departure rates must 
be the same under FIFO. This is expressed as 

r,r . .^ ,Wa(Tait))-^^, if na(t) / , , 

U^ait)) = <j ''^ uait)' ' ^ (4.19) 

if Uait) = 



w 



()4.19p uniquely determines the turning percentages at junctions with more than one outgoing 
link, and is consistent with the FIFO discipline and established route choices. It remains to 
express the path delay as the sum of finitely many link delays. If we describe path p € V as 
the following sequence of conveniently labeled arcs: 

p = {oi, 02, . . . , a-i-i, a-i, a-i+i . . . , am{p)} 

where m(p) is number of arcs in path p. 

It then follows immediately that path arrival time along p when the departure time at 
origin is t can be expressed as a composition of arc exit time functions: 

r^it) = Ta^^p^ o ... oTa^ o r^^ (t) p = {ai, a2, am(p)} GP (4.20) 
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where the operator o means decomposition, that is, / o g(x) = f(^g{x)^. 

Now the complete network loading procedure is given by ()4.14p - (j4.20p . which is interpreted 
as a well defined differential algebraic equation (DAE) system. Moreover, as well shall see in 
the next section, the (efi'ective) path delay operator defined in this way is strongly continuous 

/ \ \'P\ 

from the subset A of a Hilbert space into ( £^[to, i/] ) 

5 Existence of the DUE 

Existence results for DUE are most general if based on formulation (j3.10p . Theorem 12.51 for 
the existence of solutions of variational inequalities in topological spaces can be applied if the 
operators '^p{t, •) can be shown to be continuous and the feasible set A can be shown to be 
compact. After Section [5?T] addresses the continuity of the effective delay operator, based on 
the DNL model introduced in previously, the last obstacle to proving existence is the com- 
pactness of A, which unfortunately does not generally occur in SRDC DUE. To overcome the 
aforementioned difficulty, we will consider instead successive finite-dimensional approximations 
of A, and rely on a topological argument. Such an approach is mathematically rigorous but 
much more challenging than would be the case if A were compact in the appropriate Hilbert 
space. The topological argument and supporting infrastructure for a proof of existence are 
presented in Section 15.21 and Section 15.31 



5.1 Continuity of the effective path delay operator 

In this section, we will establish continuity of the map h i— > ^'(•, h). These results will be 
essential for the proof of existence theorem for DUE in Section 15. 3[ Notice that unlike the 
argument in Zhu and Marcotte (2000) which required a priori bound for the path flows, the 
proof provided here works for unbounded path flows and even distributions, thanks to the 
generalized Vickrey model. 

The next lemma provides the criterion for the continuity of the delay function Da{-), a £ A. 

Lemma 5.1. Consider an arc a £ A, with inflow Ua{-)- Under the generalized Vickrey model 
expressed in j^[ )- p. jg[ ), the arc delay function Da{-) is continuous if Ua{-) G >C^[io, tf]. 

Proof. Assume that Ua{-) € ^"^[to, tj], then Ua{-) G ^^[to, tf]. Therefore the cumulative 
entering vehicle count 

Ua{t) = / Ua{s)ds 



J to 

is absolutely continuous. It is straightforward to verify that the following quantity is contin- 
uous. 

qa{t) = Ua{t)-Mat- min {Ua{T)-MaT} 

T<t-Ta 

where qa{t) denotes the queue length, Ma, Ta denotes flow capacity and free flow time, re- 
spectively. By (j4.16p . the delay function Da{-) is continuous. □ 

The next lemma is a technical result that will facilitate the proof of Theorem 15.31 

Lemma 5.2. Assume /(•) : [02, 62] is continuous. Let g-ni') '■ [oi, &i] [0'2, b2], n > 1 

be a sequence of functions such that gn converges to g{-) : [ai, 61] — > [02, 62] uniformly. Then 

f{9n{-)) — > f{9{-)), n — > 00 

uniformly. 
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Proof. According to the Heine-Cantor theorem, /(•) is uniformly continuous on [a2, 62]- It 
follows that, for every e > 0, there exists 6 > such that for all yi, 1/2 S [0-2, b2], whenever 
\yi — 2/2! < the inequality 

l/(yi)-/(y2)| < e 

holds. Moreover, by the uniform convergence of g^, there exists some N > such that, for all 
n > N, we have 

\9n{x) - g{x)\ < 5, V X e [ai, 61] 

Thus, for every n > N, 

\f{gnix)) - f{gix))\ < e, VxG[ai,6i] 

□ 

Theorem 5.3. Under the network loading model described in Section [71 the effective path 

I'PI — 
delay operator ^'(t, •) : A — )• (/3^([to, tf])) , h ^ ^'(•, h) is well-defined and continuous. 

Proof. For each /i E A, the functions ^p(-, h) for sll p ^ V are uniquely determined by the 
network loading procedure. To show that the effective path delay operator is well-defined, it 

remains to show that ^'(•, h) € ^£^[to, if]^ for each /i € A. Notice that there exists an 
upper bound for the path delays regardless of the network flow profile: 



Dp{t, h) < Y.\jf Yl Q^J + Ta\, yheA,peV,te [to, tf] 

where Ma, Ta are the flow capacity and free flow travel time respectively, that are associated 
with the arc a. Recall the definition of effective path delay (j3.6p : 

^p{t, h) = Dp{t, h) + T{t + Dp{t, h) - Ta 

Since is continuous, the uniform boundedness of Dp[t, h) thus implies the uniform bound- 
edness of ^p{t, /i) for all /i E A, p G P and t G [t^, tf\. This leads to the conclusion that 

^'(•, h) € ^£^[^0, tf]^ for all /i € A. With the preceding as background, the proof of 
continuity of the effective delay operator may be given in five parts. 

Part 1. Consider first only a single link with a sequence of entry flows Uy wherein > 1 that 
converge to u in the C'^ norm. That is. 

\\uv — u\\2 = I / {uiy{t) — u{t)) dt j — > as V — > 00 

Consider the cumulative entering vehicle counts 

Uy{t) = f Uy{s)ds V > I 

7 i € [to, tf] 

U{t) = / u{s)ds 

J to 

then Uy converge to U uniformly on [to, tf]. This is due to the simple observation 

\Uy{t)-U{t)\ < I \uy{s) - u{s)\ ds < |ln,,-u||i < {to - tfY^^Wnu - u\\2 — > 

Jto 
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where || • ||i is the norm in C^[tQ, tj]. The last inequahty is a version of Jenssen's inequahty. 
Part 2. Define R{t) = U{t) — Mt, Ry{T) = U^{t) — M t. We claim the following uniform 
convergence holds 

min {i?^(T)} — y min {i?(T)} V t G [to, t/] (5.21) 

T<t T<t 

where M is the link flow capacity. Indeed, for any e > 0, by uniform convergence of Uiy, we 
can choose such that for all v > N ^ the following inequality holds 

\U^{t)-U{t)\ < e Vt G [to, tf] 

For any fixed t, if > A^, we have 

\R^{t)-R{t)\ = \U^{t)-U{t)\ < e (5.22) 

Let f = argmin^<j {R{t)}, so that, by (j5.22p . we have 

min{i?^(T)} < Ry{f) < R{f)+e = min{ii(r)} + e (5.23) 

T<t T<t 

On the other hand, let f^, = argmin^<j {R^{t)}. Then for each > A^, it must be that 

min{i?(r)} < R{f,y) < R'^'^^t^) + e = mm{RuiT)} + e (5.24) 

T<t T<t 

Taken together, ([OH]) and ([Oi]) imply 

mm{R^} - min{i?(r)} < e, M u > N 

T<t T<t 

Since t is arbitrary, the claim is demonstrated. 

Part 3. The immediate consequence of Part 2 and ()4.14p - ()4.16p is the following uniform 
convergence 

W,it) Wit), q,it) q{t), D,{t) Dit), T,{t) T{t), v ^ oo (5.25) 

for which we employ notation whose meaning is transparent. The next step is to extend such 
convergence to the whole network. Consider the sequence of departure rates hy converging to 
h in the || • ||£2 norm. By the definition (j3.5p . this implies each path fiow hp^i,{-) — > hp in the 
II • II2 norm, for all p V. A simple induction based on results established in Part 2 yields, 
as u — > 00, 

UaAi) Ua{t), WaAt) ^ Wa{t), Da,,{t) Da{t), Ta,u{t) T^it), (5.26) 

uniformly for all a (z A. 

Part 4. We will show next the uniform convergence of the path delay function Dp{-, hy) — )• 
Dp{-, h), based on (j5.26p . Recall the path exit time function (I4.20j) 

T^{t) = Ta^(p) o . . . o Taj o Ta, (t) p = {ai, a2, ■ ■ ■ , am(p)} ^ V (5.27) 

The desired uniform convergence of path delay follows from the uniform convergence T^{t) — ?> 
T^it). We start by showing that Ta2,u ° Tai,u(t) — ?> Tag o Taiit) uniformly 

For every > 1, since the inflow of arc a2 is square-integrable, Ta2,u{') is continuous by 
Lemma l5.ll This means that Ta2{-) is also continuous since it is the uniform limit of r^j, ;/(•)• 
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Lemma E2] then implies that (Taj^jy(-)) converges uniformly to (''"ai 

(•)), that is, for any 

e > 0, there exists an A^i > such that for all v > Ni, 

{t))\ < e/2, yt£ [to, tf] 
Moreover, there exists some > such that for all ly > N2, 

y 0.2,1' 

{t)-Ta,{t)\ < e/2, yt£ [to,tf] 

Now let A'o = max{A^i, For any v > n, and all t G [to, tf], 

it)) m 

it)) (t))\ + \Ta2{Ta^At)) m 

< e/2 + e/2 = e 

This shows the desired uniform convergence Ta2,u ° Tai,u{t) o Ta^{t). 

The uniform convergence t^(-) — >■ T^i ) follows immediately by (|5.27p and induction. As a 
result, we obtain the uniform convergence of path delay 

Dp{-, K) — > Dp{-, h), V — > 00 

Part 5. Finally, recall the definition of the effective delay 

^(t, h) = Dp{t, h) + F{t + Dp{t, h) - Ta) 

Note that J-{-) is continuous, the following uniform convergence follows by Lemma 15.21 

T{t + Dp{t, K) - Ta) T[t + Dp{t, h)-TA), u ^ 00 

We conclude that the effective delay ^p{-, hy) converges uniformly to ^p{-, h). The desired 
convergence in the || • ||£2 norm now follows since the interval [t^, tf] is compact. □ 

5.2 Alternative definition of effective path delay 

The integrals employed in defining the feasible domain (|3.8p are not enough to assure bounded 
path flows hp, p ^V. This observation is the fundamental hurdle to providing existence of the 
DUE solution. One of the main accomplishments of this paper is to address the boundedness of 
path flows not only for the proof of existence result but also for future analysis and estimation 
of network flows. In this section, we will present an alternative formulation of the effective 
path delay ^p{t, h), where that alternative formulation will facilitate our analysis leading to 
the proof of our main result. Theorem 15. 6i 
Recall the effective delay operator 

^p{t, h) = Dp{t, h) + T{t + Dp{t, h) - Ta) (5.28) 

In order to simplify our analysis, it is convenient to rewrite (I5.28P in a slightly different form. 
In particular, for each 0-D pair (i, j) € W, let us introduce the cost functions : [to, tf] — )• 

and ipij{-) '■ [to, tf] which measure the travel costs at the origin and destination of 

the path, respectively. More precisely, ^ij(-) is a function of departure time, while V'ij(') is a 
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function of arrival time. Fix any vector of path flows /i G A, recall the path exit time function 
T^it) = t + Dp{t, h). Then (j5.28p can be equivalently written as 

^pit,h) = -t + TP{t)+T{TP{t)-TA) 

= (j)ij{t) + 'ipij[Tp{t, h)) 

where 

(t>ij{t) = -t, i,,j{t) = t + T{t-TA) (5.29) 

Remark 5.4. In Bressan and Han (2011, 2012), the travel costs are measured in terms of 
4'ij{')j i^iji'); iiT- other words, the general effective delay 115. 28\) can he alternatively evaluated 
as a sum of costs at the beginning and end of journey of each driver. 

In Section 15. 3( we will exploit the alternative representative of effective path delay (cost) to 
establish existence. 

To prepare for the existence proof, consider a general network G{N, A). Associate with 
each 0-D pair (z, j) G W the pair of cost functions 0ij(-) and V'ij(') as defined in (I5.29p . We 
make the following two assumptions on (pij{-), ipiji') ^^^^ the underlying traffic flow model. 

Al. For all, ipij{-) is continuously differentiable on [Iq, tf]. In addition, 

j^Ajit) > 0, t G [to, tf] (5.30) 



A2. Each link a G A in the network has a flow capacity /^"^ < oo 

Note that Al stipulates that each arrival cost is strictly increasing, while A2 applies 

to all traffic flow models that assume a flow capacity for each network link. Notice that by 
()5.29p . Al amounts to requiring that 

Al'. J^{-) is continuously differentiable and satisfies J''{s) > —1 on [tg, tf]. 

Remark 5.5. In view of the preceding assumptions, we are prompted to define the following: 



d_ 

{i,f)<iwte[t(i,tj\ dt 



i^'min = , , H'^ij^^) > ^ (5.31) 



pmax ^ max/^"^ < +C50 (5.32) 

Note that V^mm > ^ follows from i5.3U\) and the fact that ^V'ij(') continuous on the compact 
interval [to, tf]. 

5.3 Existence of solution to the variational inequality 

The classical result explained by Theorem 12.51 will be the key ingredient for the proof of 
existence of the DUE solution. Using the same notation as in Theorem 12.51 the underlying 
topological vector space E will be instantiated by (>C^([to, tf])) , which is a locally convex 

I -Tjl 

topological vector space. The dual space E* will be again (>C^([to, tf])) 
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Theorem 5.6. (Existence of DUE) Let assumption (Al'), (A2) hold. In addition, assume 
that the effective delay operator ^' : A — )• (£^[to, i/]) is continuous. Then the dynamic user 
equilibrium problem as in Definition \3.1\ has a solution. 

Proof. The proof is divided into four parts. 

Part 1. Our strategy for demonstrating existence is to adapt Theorem 12.51 to the locally 
convex topological vector space (£^([to, tf])) , and its subset A. By assumption, the map 

h 1-^ ^'(•, h) is continuous from A to the space of (/3^([to, tf])) ■ If ^ were compact and 
convex, we would immediately demonstrate the desired result. However, A is bounded, closed 
and convex, but not compact in (£^([to, tf])) 

Part 2. We will instead employ finite-dimensional approximations of A. In order to 
proceed, consider for each n > 1 the uniform partition of interval [to, tj] with 2" sub-intervals 

to = t'^ < t^ < t'^ ... < t^" = tf 
i_ i_i ^ t_^-to i = i 2" 

Then consider the following sequence of finite-dimensional subsets 

An = {h£A: /ip(-) is constant on f), peV] C A (5.33) 

I'P I 

We claim that for each n > 1, A„ is compact and convex in {C\[to, tf])y' Indeed, given 
any h"''^, /i"'^ € A„, and a € [0, 1], ah"'^ -|- (1 — a) /i"'^ is clearly nonnegative and constant 
on each f] fori = 1, . . . , 2". In addition, for any origin-destination pair (i, j) € W, by 

definition (j3.8|) . 



pen, 



This verifies that A„ is convex. To see compactness, we define the map /i : 

(«i,p : 1 < i < 2", p gV) where each vector {ai^p, . . . , a2",p) is the coordinate of hp under the 

natural basis {e^}?^]^, where 



eiit) 



t e [f-i, f 

else 



Clear, the map fi is one-to-one. Now consider any sequence j/>i C An, and the sequence 

of their images {^(/i"^''^) }^>i C 3fJ+ By {M(/i"^''^}^>i is uniformly bounded by the 

following quantity 

2 Qij 



max ■ 



i,j tf — to 

By the Bolzano- Weierstrass theorem, there exists a convergent subsequence {//(/i"^*^'-*)}^,^-^. 
By construction, the subsequence {h"'^'^'^} ^i^i must converge uniformly to some /i". In view 
of the compact interval [to, tf], we conclude that this convergence is also in the || • ||£2 norm. 
It now remains to show /i" € A„, then our claim follows from sequential compactness of A„. 
Clearly /i" > and is constant on the sub-intervals t*), i = 1, . . . ,2". Moreover, since 

^71(1/') uniformly bounded, by the Dominated Convergence Theorem, 

Qij = lim y f\<'''\t)dt = y ['^h;{t)dt 
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This implies /i" G A„. Thereby, the claim is substantiated. 

Part 3. For each n > 1, apply Theorem 12.51 to A„ and obtain h^'* € A„ such that 

^^'(•, /i"'*), > 0, V/i"gA„ (5.34) 

Where the ( , ) is the inner product (duality) defined in ()3.4p . It is easy to observe that ()5.34p 
implies that if hp'*{t) > 0, t G [P , then 

/ *(t, /i"'*) = min / ^(t,h''^*)dt (5.35) 

Jti 0<k<2"Jtk 

Recall ()5.3ip . ()5.32p . choose any constant M such that 

M > (5.36) 



mm 

71, * / 



We then claim that for all n > 1, there must hold hp [t) < A4 for ail t G [Iq, tf], p & V. 
Otherwise, assume there exists some p £ V, some u > 1 and some < i <2'^ with 

/i^'*(t) = fi > M te[f, f+\ 

By choosing [Lq, tf] large enough, we can assume that i > 1. Now consider the interval 
f], and the quantity sup ^{t, H"'*), by possibly modifying the value of the function 

^(■, /i"'*) at one point, we can obtain t* G t*] such that 

^{t*, /i^'*) = sup ^{t, /i^'*) 

Now, let Tp{t, K"'*) = t -\- Dp{t, K"'*) be the arrival time function. According to FIFO and 
assumption (A2), we deduce that for all t G [t*, t*"*"^], 

rp{t,h^n-rp{t\h^n > ^^j^ (5.37) 
this implies, together with (I5.3ip that 

4^^,{Tp{t,h^n)-Mrp{t.h^n) > 4'Ln{rv{t.^"n-rr>{t\li"n) > ^'min-^^p^ (5.38) 
dOSll and (fOOl) imply 

^p{t, h^n - %it\ h^n > yt^it - n -{t- n, v t g [f , ^ +i] (5.39) 

Integrating (I5.39P from f to t*"^^ and a simple calculation yield 

(fi+l_fi\2 11 (fi+^ _ fi\ f _|_ \ 

/ *(t, /i'''*)dt-(f+i-f)^(t*, /i^'*) > ^- i^.^wf^i^ li.ff- ^ 1 

(5.40) 

Notice since t* G f ], we deduce i* - {f + t'+^)/2 > -3/2(t*+i - f). Km then becomes 
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This yields the following contradiction to (j5.35p and hence (j5.34p 

/ ^{t, h'''*)dt > / ^pit, h'''*)dt 

Part 4. By previous steps, we have obtained a uniformly bounded sequence of vector- 
valued functions {^5^}„>^ satisfying (j5.34p . By taking a subsequence, we can assume the weak 
convergence in both £^-space and £^-space 

/i"'* — > h*, n oo 

for some h* £ A. We claim that for all h £ A, 

Indeed, given any h £ A, there exists piecewise-constant approximation {/i„ G A„}„>i, which 
converges to h both point-wise and in the || • ||£2 norm. According to (I5.34p . we have 

/i"'*), > (5.42) 

Notice that the map h i— )• ^p{-, h) is continuous with respect to the || • ||£2 norm, by the 
continuity of the inner product, we pass (j5.42p to the limit 

□ 

The next theorem establishes the existence of SRDC DUE with generalized Vickrey model, 
which is an immediate consequence of Theorem 15.31 and Theorem 15.61 

Theorem 5.7. Let assumption (Al ') hold, then the dynamic user equilibrium as in Definition 
\3.1\ with generalized Vickrey model has a solution. 

Remark 5.8. Notice that the assumption A2 is satisfied by the Vickrey 's model and hence 
the GVM. Thus in the statement of the theorem, A 2 is omitted. 



6 Conclusion 

We have established the existence of the continuous-time simultaneous route-and-departure 
choice DUE for the generalized Vickrey model (GVM), a generalization of Vickrey's model, 
and plausible regularity conditions that are easy to check and rather weak. It is significant 
that ours is the first DUE existence result without the a priori bounding of departure rates 
(path fiows) and the most general constraint relating path flows to the trip table. In fact, our 
method of proof successfully overcomes two major hurdles that have stymied other researchers: 

1. the set of feasible flows A is intrinsically non-compact in >C^-space as well as £^-space; 
and 

2. a direct topological argument requires a priori bounds for the path flows, where those 
bounds do not arise from any behavioral argument or theory. 

Theorem 15.61 is a general result that subsumes all SRDC DUE models regardless of the arc 
dynamics, flow propagation and arc delay function employed, as long as the effective delay 
operator is continuous. 
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